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Objectives 


The  main  challenge  in  structural  health  monitoring  is  to  identify  the  location  and  severity  of 
small-sized  damage  at  an  early  stage.  The  two  classical  and  popular  approaches,  the  vibration- 
based  methods  and  the  ultrasonic  wave-based  methods,  both  have  limitations.  This  research 
explores  damage  identification  via  advancing  a  third  type  of  approach;  high-frequency  harmonic 
excitation-based  self-sensing,  i.e.,  the  impedance-based  approach.  Specifically,  we  propose  to 
create  a  new  concept  of  adaptive  high-frequency  piezoelectric  self-sensing  interrogation  by 
means  of  tunable  circuitry  integration  to  the  piezoelectric  transducer.  The  underlying  hypothesis 
is  that,  by  tuning  online  the  circuitry  elements  properly,  both  the  quality  and  the  quantity  of  high- 
frequency  admittance  (i.e.,  the  inverse  of  impedance)  measurement  data  can  be  greatly  increased 
for  damage  detection  purpose,  which  provides  a  foundation  for  the  subsequent  development  of 
new  damage  identification  algorithms.  The  overarching  objective  is  to  fundamentally  enhance 
the  identification  accuracy  and  confidence-level  by  the  concurrent  advancement  of  sensing 
hardware  and  inverse  identification  algorithms. 

Summary  of  Efforts 

In  this  project,  a  series  of  efforts  are  undertaken  to  validate  the  original  hypothesis  and  to 
formulate  a  systematic  methodology  of  high-frequency  harmonic  response-based  damage 
identification,  which  includes: 

•  At  the  sensing  mechanism  level,  new  concept  of  enhanced  piezoelectric  admittance  sensing 
using  circuitry  integration  is  developed,  and  sensor  design  guidelines  are  provided. 

•  At  the  inverse  identification  methodology  level,  new  algorithms  of  damage  identification  are 
synthesized,  which  can  take  full  advantage  of  the  new  sensing  mechanism. 

These  efforts  have  thoroughly  addressed  the  research  issues  outlined  in  the  project  proposal,  and 
yielded  a  complete  methodology  of  adaptive  high-frequency  piezoelectric  self-sensing 
interrogation. 


Descriptions  of  Accomplishments  and  New  Findings 

The  timely  detection  of  structural  faults  in  aircraft  components  has  obvious  significance  to 
the  Air  Force.  The  past  decade  has  witnessed  ever-increasing  enthusiasm  on  developing 
structural  health  monitoring  (SHM)  systems  that  can  provide  near  real-time  detection  of 
structural  damage.  Currently  there  are  two  widely  studied  classes  of  methods  in  SHM  utilizing 
structural  dynamic  responses:  the  vibration-based  and  the  wave  propagation-based  methods.  In 
both  classes  of  methods,  embedded  piezoelectric  transducers  have  been  widely  used  to  enable 
online  monitoring  of  structural  health.  The  piezoelectric  transducers  possess  two-way  electro¬ 
mechanical  coupling  (i.e.,  deforming  when  subject  to  electrical  excitation,  and  generating 
voltage/current  when  subject  to  deformation)  and  behave  electrically  as  capacitors.  Owing  to 
these  features,  when  we  integrate  (bond/embed)  a  piezoelectric  transducer  to  a  structure,  the 
electrical  impedance/admittance  of  the  transducer  is  directly  related  to  the  impedance  of  the 
underlying  structure.  This  has  led  to  the  recent  interest  in  developing  another  class  of  embedded 
SHM  methods,  the  piezoelectric  impedance/admittance-based  methods. 
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The  piezoelectric  impedance-  or  admittance-based  approach  appears  to  be  quite  promising 
because  of  the  following  reasons;  1)  Compared  to  the  vibration-based  methods,  the  piezoelectric 
impedance/admittance  can  be  extracted  in  high-frequency  range  (e.g.,  higher  than  30  kHz), 
which  is  important  for  detecting  small-sized  damage;  and  2)  Compared  to  the  wave-propagation 
methods,  the  impedance/admittance  is  extracted  from  the  stationary  responses  (i.e.,  harmonic 
responses),  and  thus  is  better  suited  for  damage  identification  (i.e.,  identifying  both  the  damage 
location  and  severity).  On  the  other  hand,  it  is  worth  noting  that  while  the  piezoelectric 
impedance  or  admittance-based  approach  holds  great  potential,  the  reality  is  that  such  potential 
has  not  been  fully  realized.  There  appear  to  be  severe  bottlenecks:  1)  While  actual  detection 
sensitivity  hinges  upon  the  scale  of  difference  in  the  measurements  upon  the  occurrence  of 
damage,  the  response  anomaly  can  still  be  buried  in  the  measurement  noise  because  different 
damage  profiles  affect  the  responses  differently;  2)  These  methods  generally  have  limitation  of 
inaccuracy  due  to  under-determined  inverse  problem  for  damage  prediction,  especially  under 
noise  influences;  and  3)  Current  practices  mainly  resort  to  ad  hoc  damage  metrics,  i.e.,  to  declare 
damage  occurrence  if  the  difference  between  the  online  impedance/admittance  measurement  and 
the  healthy  baseline  exceeds  a  certain  threshold,  and  are  still  unable  to  accurately  identify 
damage  location  and  severity. 

Based  upon  the  state-of-the-art  review  and  the  Pis’  experience,  we  have  recognized  that  some 
characteristics  of  the  piezoelectric  circuitry  can  be  utilized  to  revolutionize  the  piezoelectric 
impedance/admittance  sensing.  In  particular,  we  propose  to  integrate  tunable  circuitry  with  the 
piezoelectric  transducer  to  realize  dynamical  tailoring  of  the  sensing  system  to  fundamentally 
enhance  the  sensitivity  and  robustness.  The  major  accomplishments  and  new  findings  garnered 
throughout  this  project  are  summarized  in  the  following  sections. 

Part  I.  Integration  of  adaptive  circuitry  to  enhance  the  quality  and  quantity  of  impedance 
measurements 

In  this  first  part  of  research,  we  explore  the  fundamental  sensing  mechanism  enhancement  by 
circuitry  integration.  For  simplicity  in  discussion,  here  we  use  a  one  degree-of-freedom  (DOF) 
model  to  represent  the  interested  mode  of  the  mechanical  structure.  The  dynamic  equations  of 
the  structure  coupled  with  the  piezoelectric  transducer  without  and  with  the  inductance  can  be 
derived  as 


mq  +  cq  +  kq  +  k^Q-0,  RQ  +  k2Q  +  k^q  =  V-  (without  inductance)  (la,b) 

mij  +  cq  +  kq  +  k^Q-O ,  LQ  + RQ  + k2Q  + k^q -V.  (with  inductance)  (2a,b) 

where  q  is  the  mechanical  displacement,  Q  is  the  electrical  charge  in  the  circuit,  m,  c,  and  k  are 
the  mass,  damping,  and  stiffness,  respectively;  k^^HC^  is  the  inverse  of  the  piezoelectric 

capacitance,  and  FJ  is  the  excitation  voltage.  Both  systems  include  a  resistor  that  enables  the 
direct  measurement  of  current  resulted  from  the  voltage  excitation.  The  piezoelectric  admittances 
can  be  expressed  as,  respectively. 
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(without  inductance) 


(3) 


/  icoi^-mco^  +icco  +  k^ 

’’  V.  {iRco  +  k2)[-mo/  +  icco  + 

I  io)i-mo)^+icG>  +  k\ 

7  =  —  =  - - - - - r - -  (with  inductance)  (4) 

Vi  [-Lor  +  iRco  +  k^  \  [-mor  +  ico)  +  k\-  A:, 

While  the  impedance  and  the  admittance  are  inverse  with  each  other,  in  this  research  we  focus  on 
the  admittanee,  beeause,  as  shown  in  the  comparison  between  Eqs.  (3)  and  (4),  the  admittanee 
measurements  ean  be  greatly  impacted  by  the  eircuitry  dynamies.  Indeed,  one  ean  see  that  Eq.  (4) 
is  a  fourth-order  system  with  two  resonanees,  while  Eq.  (3)  is  essentially  a  seeond-order  system 
with  one  resonance  sinee,  compared  to  k^,  Rco  is  very  small.  Obviously,  one  major  difference 
between  the  new  sensor  with  the  induetance  (Eq.  (4))  and  the  traditional  approaeh  (Eq.  (3))  is 
that  the  induetive  circuitry  introduces  an  additional  DOE,  which  leads  to  a  new  resonant  effect  in 
the  admittanee  curve.  It  is  well  known  that  the  damage-induced  ehanges  are  most  significant 
around  the  resonant  peaks.  Our  preliminary  investigation  indieated: 

•  The  circuitry  resonance  effect,  if  properly  utilized,  may  amplify,  under  the  same  level  of 
voltage  exeitation,  the  magnitude  of  the  admittance.  Eurthermore,  it  may  even  amplify  the 
damage-induced  admittance  change. 

•  Through  tuning  the  inductanee  to  different  values,  we  ean  obtain  a  mueh  enlarged  dataset  for 
damage  deteetion,  which  consists  of  a  family  of  admittance  curves  (under  different 
inductanee  values)  as  eompared  to  the  original  single  admittance  eurve  obtained  by  the  fixed 
sensor  without  circuit. 

The  new  idea  introduced  in  this  part  of  researeh  is  to  integrate  a  negative  capaeitance  into  the 
cireuit.  In  this  ease,  we  can  obtain  the  eleetrical  admittance  of  the  enhanced  series  circuitry, 
normalized  by  the  eleetrical  admittance  of  the  stand-alone  piezoelectric  transducer  (before  it  is 
integrated  onto  the  meehanical  structure),  whieh  is  Yp  =  icolk^ ,  as 


kj  [-mof  +  ica>  +  k) 

Yp  i^-Lof  +  iRo)  +  ^2  +  ico)  +  k)  -  k^ 


The  above  expression  can  be  can  be  written  as 

i>  5^[\  +  2giLl-Q}) 

7p  ~  (l  +  IgiQ.  -  Q" )( +  rSHa  -  Q" )  - 


(5) 


(6) 


where 


co„ 


c 

2ma)^ 

m 


K 


(7a-g) 
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Here  ,  co^,  Q. ,  5 ,  g ,  r  and  ^  are  the  meehanieal  modal  frequency,  the  electrical  circuitry 
resonant  frequency,  the  non-dimensionalized  excitation  frequency,  the  ratio  between  the 
electrical  resonant  frequency  and  the  mechanical  modal  frequency  (i.e.,  inductance  tuning  ratio), 
the  structural  damping  ratio,  the  electrical  damping/resistance  ratio,  and  the  generalized  electro¬ 
mechanical  coupling  coefficient,  respectively.  Without  loss  of  generality,  the  structural  damping 
ratio  is  set  as  ^  =  0.001 .  Clearly,  after  the  negative  capacitive  element  is  integrated,  the  inverse  of 
the  overall  capacitance  in  the  circuit  is  reduced,  and  hence  the  generalized  electro-mechanical 
coupling  coefficient  is  increased,  i.e.,  ^  . 


The  circuitry  admittance  curves  of  the  benchmark 
undamaged  beam  with  different  generalized  coupling 
coefficients  are  shown  in  Figure  1.  The  coupling  coefficients 
are  0.0387,  0.0726  and  0.1547,  which  correspond  to  the  cases 
without  the  negative  capacitance  (piezoelectric  inductive 
circuitry),  with  a  -12.39  nF  negative  capacitance,  and  with  a  - 
9.46  nF  negative  capacitance,  respectively.  Consider  the  17*^ 
mechanical  modal  frequency  (6120  Hz).  In  order  to  have  ^  =  1 
(to  achieve  large  dynamic  coupling  between  the  transducer  and 
the  host  structure),  under  these  coupling  coefficients,  the 
inductance  values  are  selected  as  76.2  mH,  21.65  mH  and  4.75 
mH  (which  are  actually  the  frequency-veering  inductance  values  under  these  three  coupling 
coefficients,  correspondingly.  For  reference,  the  admittance  curve  extracted  by  the  low-cost 
impedance  measurement  circuit  (without  the  tunable  inductance  and  negative  capacitance)  is  also 
plotted.  It  can  be  seen  that  the  larger  the  coupling  coefficient,  the  flat  the  admittance  curve  is 
whereas  the  resonant  peaks  remain  about  the  same.  For  example,  when  the  system  is  excited  at  4 
kHz,  if  we  integrate  a  negative  capacitive  element  with  a  0.0726  coupling  coefficient,  the 
admittance  magnitude  is  15.83  dB,  which  is  10.96  dB  larger  than  that  without  the  negative 
capacitance  (4.87  dB).  When  we  connect  a  negative  capacitive  element  corresponding  to  a  larger 
coupling  coefficient  =  0.1547),  the  increased  admittance  magnitude  is  22.85  dB,  compared 
with  the  result  without  the  negative  capacitance.  This  means  that  the  admittance  magnitude 
remains  at  a  high  level  in  a  much  wider  frequency  range  around  the  resonant  peaks  of  the 
admittance  curve,  when  larger  coupling  coefficient  is  selected.  Therefore,  the  negative 
capacitive  element  can  further  increase  the  signal-to-noise 
compared  with  the  preliminary  study. 

The  more  important  effect  of  the  negative  capacitive 
element  can  be  seen  in  the  change  of  admittance  signature  when 
damage  occurs.  This  is  because  the  larger  the  change  of 
admittance  under  damage  occurrence,  the  more  sensitive  the 
damage  detection  scheme  is.  In  this  research,  the  structural 
damage  is  assumed  to  cause  the  loss  of  beam  equivalent  modal 
stiffness  {k),  and  the  severity  of  damage  is  defined  as  its 
percentage  change.  After  a  1%  beam  modal  stiffness  loss 
occurs,  the  changes  of  the  admittance  under  different  coupling 


ratio  in  sensor  measurement. 


Figure  2  Circuitry  admittance 
change  curves  under  different 
generaiized  coupiing  coefficients. 


Figure  1  Circuitry  admittance 
curves  under  different  generaiized 
coupiing  coefficients. 
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coefficients  are  shown  in  Figure  2.  We  can  see  that  the  larger  the  coupling  coefficient,  the  flat 
the  curve  of  admittance  change  is  whereas  the  peaks  are  reduced  only  slightly.  For  example, 
with  the  voltage  excitation  frequency  at  5  kHz,  the  magnitude  of  admittance  change  is  increased 
21.93  dB  after  a  -12.39  nF  negative  capacitance  is  connected  (if  =  0.0726  ).  If  the  negative 
capacitance  with  a  larger  coupling  coefficient  (-9.46  nF  and  f  =  0.1 547  )  is  integrated,  the 
increased  admittance  change  at  this  frequency  is  42.77  dB.  It  is  well  known  that  in  the 
impedance/admittance  based  sensory  system,  there  exist  frequency  shifts  and  magnitude  changes 
in  the  piezoelectric  impedance  or  admittance  curves  after  the  structural  damage  occurrence,  and 
these  changes  are  most  significant  around  the  resonant  peaks.  Therefore,  the  negative  capacitive 
element  can  enrich  the  input  information  for  the  damage  detection  and  identification,  and  hence 
improve  the  detection  accuracy  and  sensitivity. 


Part  II,  Analysis  on  the  adaptive  piezoelectric  circuitry  for  data  enrichment 


Piezoelectric 

transducer 

/  ^ 


Ve  Excitation  voltage 
_ K 


This  part  of  the  research  develops  the  damage 
identification  algorithm  which  exploits  the 
enriched  admittance  information  obtained  by 
employing  the  adaptive  piezoelectric  circuitry. 

Numerical  analysis  is  performed  by  using  the 
spectral  element  method  (SEM),  and 
experimental  investigation  clearly  verified  the 
effectiveness  of  the  proposed  methodology. 

Figure  3  shows  the  schematic  of  an  example  structure  integrated  with  adaptive  piezoelectric 
circuit  for  data  enrichment.  The  piezoelectric  transducer  on  the  P^th  element  of  the  beam 
structure  is  used  for  piezoelectric  impedance  measurement,  and  another  transducer  with  tunable 
inductor  is  implemented  on  the  P2th  element  for  data  enrichment. 


Figure  3  Illustration  of  the  structure  integrated  with 
the  adaptive  piezoelectric  circuitry. 


First,  the  transducer  on  the  P^th  element  is  only  considered  in  the  analysis.  As  we  assume  that 
the  structural  damage  is  represented  by  the  element  bending  stiffness  reduction,  the  piezoelectric 
admittance  of  the  healthy  and  damaged  structures  can  be  formulated  as  T(a)),  T(a))  by  using 
spectral  element  model,  respectively. 


r(m)  =  +  1] 

+S-Hco)[Ef=iCfK,iCidjS-i(co)]<I>^  +  1] 


4s 


iCx)RCn 


(8) 

(9) 

(10a,b) 


where  £^3,  Ip,  bp,  hp,  Ep,  and  Cp  are  the  permittivity,  length,  width,  thickness,  piezoelectric 
coefficient,  elastic  modulus,  and  capacitance  of  the  piezoelectric  transducer,  respectively,  hp  is 
the  thickness  of  the  beam,  o)  is  the  excitation  frequency,  and  R  is  the  resistance.  S(a))  is  the 
dynamic  stiffness  matrix  of  the  beam  structure  formulated  by  the  SEM.  Op  indicates  the  location 
of  the  transducer  at  the  P^th  element,  and  dj  is  the  stiffness  reduction  of  the  ith  element.  Cj  is  the 
connectivity  matrix  for  assembling  elemental  matrices,  and  Kgj  is  the  ith  elemental  stiffness 
matrix  formulated  by  finite  element  method,  which  is  used  to  approximate  the  partial  derivative 
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of  the  damaged  stiffness  matrix  inverse  S“^(a)).  The  first-order  sensitivity  equation  that  direetly 
relates  the  differenee  of  the  undamaged  and  damaged  piezoeleetrie  admittanee  AK  with  damage 
indiees  ean  be  derived  as  the  following. 


AY^^K^K2%[S-\aj)[Y.UcjK,iCidi]S-Kaj)]<p; 

=  Gxd  (11) 

AY  =  [f  (mi),  YM . YMf  -  [T(mi),  YM . Y^V  (12a) 

Gji  =  ^K^K2<Pp[S-\a)j)CjK,tCtS-\a)j)]<Pl  i=l,  2,  N  andj=l,  2,  m  (12b) 

d  =  [di,  d.2, ... ,  dj^  (12e) 


where  m  is  the  number  of  admittanee  measurements  and  N  is  the  number  of  elements  in  the 
model.  The  loeation  and  severity  of  damage  ean  be  identified  by  solving  the  damage  index  veetor 
d.  However,  in  most  praetieal  oases,  d  oannot  be  uniquely  obtained  and  beoomes  highly  sensitive 
to  measurement  errors  sinoe  the  sensitivity  matrix  G  is  rank  defioient,  i.e.  m  <  N . 


This  oonoern  of  under-determined  inverse  problem  is  addressed  by  implementing  the  oonoept 
of  adaptive  piezoeleetrie  oirouitry  for  impedanoe  data  enriohment.  As  we  implement  additional 
piezoeleetrie  transduoer  with  tunable  induotor  L  as  shown  in  Figure  3,  a  new  degree  of  freedom 
is  oreated  by  the  oirouit  elements  into  the  integrated  system.  By  seleotively  tuning  the  induotanoe 
values,  the  dynamics  of  the  coupled  system  can  be  continuously  and  systematically  modified, 
which  enables  to  obtain  different  impedance  responses  with  respect  to  same  structural  damage. 
The  modified  stiffness  matrix  of  the  integrated  system  can  be  obtained  as  following. 


K 


lA  — 


s  =  s  — 

*^new 

)  (ud'LCp.^') 
LC-p2) 


(13) 

(14) 


where  Op^and  Cp^  indicate  the  location  and  capacitance  of  the  additional  piezoelectric  transducer, 
respectively.  The  inverse  equation  for  the  integrated  system  can  be  similarly  formulated. 

~  Gnew  ^  d  (15) 

AFjjgyi,((u)  l^new(^)  l^Tiew(^)  (16a) 

Gnewji  =  ^K-i^K2<Pp[S-^y,{a)j)CjKgiCiS-^y,{a)j)]<P^  i=l,  2,  ...,N  andj=l,  2,  ...,  m  (16b) 


When  the  inductance  is  tuned  to  form  a  sequence  from  to  Lp,  we  can  derive  a  family  of  P 
different  inverse  equations  and  augment  them  in  the  matrix  form  as  following. 


Slim  ~  Gsum  ^  d 


(17) 
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AY  =  \ 

sum  ] 

I  : 

^ '  ^sum 

^newC  ^l) 
^new(  ^2) 

(18a,b) 

I 

^AKjjgyj,(  Tp)J 

-^newC  l-‘P^  - 

Equation  (17)  shows  the  merit  of  using  adaptive  piezoeleetrie  eireuitry.  The  original  sensitivity 
matrix  {m  by  N)  is  rank  defieient,  i.e.  m  <  N.  As  P  different  induetanee  values  are  applied 
to  the  eireuitry,  the  number  of  impedanee  measurements  and  the  sensitivity  matrix  is  greatly 
enriehed  up  to  P  x  m,  thus  the  originally  ill-posed  inverse  problem  is  signifieantly  improved. 

Numerical  analysis  on  the  damaee  identification  usins  adaptive  piezoelectric  circuitry 

In  this  seetion,  damage  identifieation  studies  are  provided  to  examine  the  viability  of  the 
proposed  approaeh  for  impedanee  data  enriehment.  The  eonfiguration  of  the  illustrative  system  is 
given  in  Figure  3.  A  fixed- fixed  beam  is  divided  into  31  elements,  and  piezoeleetrie  transdueers 
are  integrated  at  the  3'^‘*  and  2E*  element.  The  dynamies  of  the  integrated  system  is  altered  by  the 
adaptive  piezoeleetrie  eireuitry  at  the  3'^‘^  element,  and  the  other  transdueer  is  utilized  for 
impedanee  measurement.  Some  relevant  system  parameters  are  presented  in  Table  1.  In  this 
analysis,  Tikhonov  regularization  is  employed  to  solve  the  inverse  problems. 

In  the  first  ease  study,  the  damage  is  assumed  as  10%  stiffness  reduetion  on  the  13th  element. 
Random  noise  of  62  dB  signal-to-noise  ratio  (SNR)  is  added  to  the  impedanee  response.  Sinee 
the  damage  effeet  is  more  signifieant  near  resonanee  peaks,  the  impedanee  ehanges  around  eaeh 
resonanee  peak  are  utilized  for  the  damage  identifieation.  Then,  the  induetanees  are  tuned  onee 
for  eaeh  resonanee  peak  to  enrieh  the  impedanee  measurements.  Table  2  provides  8  induetanees 
seleeted  for  all  resonant  frequeneies  in  the  frequeney  range  of  5  kHz  to  1 1  kHz. 


Table  1  System  parameters 


Beam  structure 

Piezoelectric  material 

Length  X  Thickness  X  Width,  mm 

607.8  X  7.62  X  3.175 

Length  X  Thickness  X  Width,  mm 

17  X  7.62  X  0.191 

Young’s  modulus,  GPa 

73.4 

Young’s  modulus,  GPa 

66 

Density,  kg/m^ 

2780 

Density,  kg/m^ 

7800 

Loss  factor,  % 

0.15 

Permittivity,  £33  ,  F/m 

1.6X  10“® 

Poisson’s  ratio 

0.33 

Piezoelectric  constant,  d3i,  mfW 

-190X  10-12 

Table  2  Inductance  values  selected  for  each  resonance  peak 


Resonant  frequency,  kHz 

5.39 

6.05 

6.75 

7.49 

8.26 

9.06 

9.91 

10.79 

Inductance,  mH 

134.6 

106.4 

85.1 

68.9 

56.4 

46.6 

38.8 

32.6 

As  we  utilize  these  enriehed  impedanee  measurements,  the  total  number  of  piezoeleetrie 
impedanee  data  set  ean  be  inereased  by  8  times  in  this  ease  study.  A  eomparison  of  the  damage 
predietion  results  with  and  without  this  impedanee  data  enriehment  is  provided  in  Figure  4.  The 
horizontal  axis  is  the  element  number  and  the  damage  index,  stiffness  reduetion  of  the  element, 
is  provided  along  the  vertieal  axis.  The  dotted  line  shows  the  aetual  damage  distribution.  It  ean 
be  observed  that  the  result  obtained  by  traditional  method  without  data  enriehment  eannot 
aeeurately  identify  damage  eharaeteristies  due  to  its  severely  ill-posed  inverse  equation,  even 
with  a  low  level  of  additive  noise  (62  dB  SNR).  However,  when  the  proposed  method  of 
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adaptive  circuitry  is  applied,  the 
damage  identification  result  indicates 
approximately  10%  stiffness  reduction 
at  the  13*  element.  The  root  mean 
square  deviation  (RMSD)  between  the 
predicted  stiffness  loss  d  and  the  actual 
damage  d°  is  employed  to  quantify  the 


(t.) 
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prediction  error  and  compare  the  Figure  4  Identification  of  damage  on  the  13'"  element  of  the 
performance  of  damage  identification.  example  beam  structure  with  10%  stiffness  reduction. 


RMSD  =  /m^df  (19) 

where  di^di  are  the  predicted  and  actual  damage  indices  of  the  zth  element,  respectively.  When 
the  adaptive  circuitry  is  applied,  the  RMSD  is  reduced  from  72  %  to  29.2  %. 


Another  numerical  case  study  is  „  o#rr 


!  th  o  -4 
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I  Without  Data  Enrichment 

RMSD:  85.6% 
I  With  Data  Enrichment 
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conducted  with  multiple  damages.  The  “ 
structural  damages  are  assumed  to  be  5  %| 
and  8  %  stiffness  reductions  at  the  13 
and  18*  elements,  respectively.  In  this  , 
case,  higher  level  of  noise  (36  dB  SNR)  ' 
is  added  in  the  impedance  response. 

Figure  5  shows  the  damage  prediction 
result.  Without  the  adaptive  circuitry,  the 
damaged  elements  are  virtually 

indistinguishable,  and  the  result  falsely  predicts  the  largest  two  damages  on  the  15*  and  18 
elements  with  approximately  2  %  stiffness  reduction.  Flowever,  as  inductance  tunings  are  applied, 
both  damage  indices  at  the  13*  and  18*  elements  clearly  indicate  the  simulated  damages  in  the 
structure.  The  damage  identified  by  using  the  traditional  method  has  85.6  %  of  the  RMSD  error, 
whereas  the  proposed  approach  shows  the  RMSD  of  54.7  %. 


Figure  5  Multiple  damage  identification  on  the  13  , 18 
elements  of  the  beam  with  5%,  8%  stiffness  reduction, 
respectively. 


th 


Investisation  on  the  number  of  inductance  tuninss 

Singular  value  decomposition  is  conducted  to  examine  how  the  number  of  inductance  tunings 
affects  the  condition  of  the  augmented  sensitivity  matrix.  A  case  study  is  performed  at  the 
frequency  range  of  5  kHz  to  11  kHz,  which  contains  8  resonance  frequencies,  and  here  we 
increased  the  number  of  inductance  tunings  to  20.  The  inductance  values  are  first  selected  to 
match  the  electrical  resonance  of  the  adaptive  piezoelectric  circuitry  with  each  structural 
resonance  of  the  host  structure  in  the  interested  frequency  range  as  shown  in  Table  1.  The  next 
tuning  values  are  deviated  ±  1  %  from  the  previously  selected  inductances  for  resonances,  and 
then  ±  2  %,  and  so  on. 


Figure  6(a)  presents  the  logarithm  of  the  singular  values  of  the  augmented  sensitivity  matrix 
^sum  in  equation  (17).  As  shown  in  the  figure,  the  singular  values  are  greatly  increased  overall 
when  the  first  8  inductances  are  applied  for  impedance  data  enrichment,  and  they  increased  even 
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more  as  the  number  of  inductance  tunings  is  increased  to  20.  This  implies  that  the  ill-posed 
sensitivity  matrix  can  be  improved  as  more  inductance  tunings  are  applied.  As  a  result,  the 
damage  prediction  error  (RMSD)  decreases  when  more  inductance  tunings  are  introduced  in  the 
system  as  shown  in  Figure  6(b).  We  can  notice  that  the  damage  identification  error  dramatically 
decreases  for  the  first  several  inductance  tunings,  and  then  decreases  gradually  afterwards. 

(a)  2 
0 

2-2 

°  -4 
-6 
-8 

5  10  15  20  25  30  5  10  15  20 

SVD  Number  Number  of  Inductance  Tunings 

Figure  6  (a)  Singular  values  of  the  augmented  sensitivity  matrix  when  different  numbers  of 
inductance  tunings  are  employed,  (b)  Damage  prediction  error  (RMSD)  with  respect  to  the  number 
of  inductance  tunings. 

Same  analysis  is  conducted  on  other  cases  with  4  resonances  in  the  frequency  range  from  8 
kHz  to  11  kHz,  and  12  resonances  in  the  range  of  4  kHz  to  13  kHz  to  investigate  whether  the 
damage  prediction  errors  decrease  similarly  in  other  frequency  ranges.  All  three  case  studies  are 
repeated  5  times,  and  the  average  damage  prediction  errors  for  each  case  are  compared  in  Figure 
7.  The  highlighted  boxes  in  the  figure  show  that  the  RMSD  starts  to  decrease  gradually  from 
near  the  4*,  8*,  and  12‘^  inductance  tuning  augmentations  when  there  are  4,  8,  and  12  resonances 
in  the  interested  frequency  ranges,  respectively.  These  results  indicate  that  there  appears  to  be  a 
relation  between  the  number  of  inductance  tunings  and  resonances  employed  for  damage 
identification.  This  can  be  explained  by  examining  how  the  inductance  values  are  selected,  and 
how  they  affect  the  inverse  equation  from  the  case 
study  of  8  resonances  as  an  example.  As  described 
previously,  the  first  8  inductance  values  are  tuned  to 
each  structural  resonance,  and  then  the  next 
inductance  values  are  varied  from  these  values. 

When  the  inductances  are  tuned  for  each  resonant 
frequency,  the  piezoelectric  impedances  of  the 
integrated  system  are  altered  distinctively  from  each 
other.  As  these  responses  are  arranged  in  the  form  of 
equation  (17),  we  can  fully  exploit  the  merit  of  data 
enrichment  by  augmenting  unique  and  independent 
information  into  the  sensitivity  matrix.  Therefore,  we 
can  greatly  improve  the  ill-posedness  of  the  inverse 
equation,  which  results  in  notable  decrease  of  the 
damage  prediction  error.  However,  the  next  groups 
of  inductances  are  variations  from  the  first  group  of 
8  inductances,  thus  the  corresponding  impedance 
responses  are  variations  of  the  first  group  of 
impedance  responses.  This  set  of  inductances  adds 
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Figure  7  Damage  prediction  error  (RMSD) 
comparisons  for  frequency  ranges  with 
different  number  of  resonances. 
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analogous  information  to  the  sensitivity  matrix,  and  thus  the  ill-posedness  is  improved  gradually. 
For  this  reason,  the  damage  identification  error  decreases  gradually  near  after  the  8*  inductance 
tunings.  From  these  results,  we  can  conclude  that  applying  more  inductance  tunings  can  increase 
the  accuracy  in  damage  prediction  by  improving  the  ill-posed  nature  of  the  sensitivity  matrix.  On 
the  other  hand,  it  is  suggested  that  the  inductance  can  be  tuned  once  for  each  resonance  in  the 
frequency  range  of  interest  to  efficiently  exploit  the  benefit  of  data  enrichment  in  damage 
identification. 


Experimental  Validation 

A  fixed-fixed  aluminum  beam  (Al-2024) 
discretized  into  61  elements  is  used  in  this 
analysis,  and  two  piezoelectric  transducers 
(PSI-5A4E)  are  bonded  on  the  top  surface 
of  the  29**^  and  4fi‘  elements  for  inductance 
tuning  and  impedance  measurement 
purposes,  respectively.  We  introduced 
damage  as  a  surface  notch  at  the  25**^ 
element  of  the  beam  structure  for  the  damaged  case.  The  notch  is  0.09  mm  deep  (2.8%  of  the 
beam  thickness)  and  10.4  mm  long  having  same  width  of  the  beam,  which  results  in  8.3% 
reduction  of  the  local  bending  stiffness.  The  experimental  configuration  including  the  location  of 
the  piezoelectric  transducers  and  damage  is  presented  in  Figure  8.  A  tunable  inductor  used  in  this 
study  is  realized  by  op-amps  (LM324),  resistors,  and  capacitors.  Table  3  provides  the  relevant 
parameters  and  the  seven  inductance  tunings  employed  for  data  enrichment.  The  baseline  model 
is  tailored  by  modifying  the  length,  mechanical  loss  factor  of  the  beam,  and  the  location  of  the 
piezoelectric  transducers  to  minimize  the  measurement  errors  of  the  beam  dimensions  and 
material  properties  in  modeling.  The  measurement  noise  is  42  dB  SNR,  however  considering  the 
uncertainty  in  the  baseline  model,  the  total  noise  level  can  increase.  In  this  analysis,  the  damage 
indices  are  solved  by  using  least-square  algorithms  with  negative  constraints.  The  coefficients 
are  restricted  to  be  negative  since  we  assumed  damage  as  stiffness  reduction  by  designing 
surface  notch  in  the  experiment. 
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Figure  8  Configuration  of  the  experimentai  beam 
structure  (measured  in  miiiimeters). 


Tabie  3  Parameters  of  the  tunabie  inductor,  beam  structure,  and  piezoeiectric  transducers 


Beam,  mm 

627.2  X  7.21  X  3.175 

Pi,  mm 

16.85  X  7.09  X  0.191 

P2,  mm 

16.3  X  7.09  X  0.191 

Circuitry 

element 

Ci,nF 

C2,nF 

R2,  fl 

R3,  fl 

R4,  fi 

R5,  fi 

Re,  n 

Ry,  fl 

9.788 

9.807 

4620 

4629 

4629 

4628 

4632 

4629 

L,  mH 

39.45 

40.367 

41.272 

45.32 

46.16 

63.1 

80.12 

Ri,  fi 

850 

870 

890 

980 

1000 

1371 

1750 

Figure  9  shows  the  damage  prediction  result  and  the  damage  identification  errors  with  respect 
to  the  number  of  inductance  tunings.  Without  data  enrichment,  the  result  indicates  damage  at 
false  location  (22"‘*  element)  with  4.3  %  stiffness  reduction.  Since  the  damaged  locations  may  not 
be  known  in  advance  in  practical  circumstances,  and  the  structural  health  monitoring  system  may 
not  be  able  to  determine  erroneous  results,  the  incorrect  predictions  would  mislead  the  following 
decisions  for  remedy.  However,  as  the  inductance  tunings  are  applied,  the  damage  prediction 
error  gradually  decreased  and  the  final  result  correctly  identified  the  damage  at  the  25**^  element 
with  predicted  stiffness  loss  (6.8%)  close  to  the  actual  8.3%  value.  These  results  verify  that  the 
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proposed  approach  of  data  enrichment  significantly  improves  the  damage  identification  accuracy 
under  noise  influence. 
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Figure  9  (a)  Experimental  damage  identification  on  the  25*^  element  of  the  beam  structure,  (b)  Damage 
prediction  error  (RMSD)  with  respect  to  the  number  of  inductance  tunings. 

Part  III,  Bayesian  inference-based  damage  identification  using  piezoelectric  admittance 
with  adaptive  circuitry 

An  important  element  in  the  research  of  piezoelectric  impedance  method  is  to  evaluate  the 
location  and  severity  of  damage  using  measurement  data.  While  the  spectral  element  method 
combined  with  the  adaptive  circuitry  idea  has  exhibited  excellent  performance,  for  very  complex 
structures  the  damage  identification  based  on  the  inverse  analysis  of  such  formulation  may  still 
encounter  difficulty,  because  of  high  computational  cost  and  the  usual  under-determined  nature 
of  the  problem.  That  is,  for  structures  with  complex  geometry  and  very  large  size,  the  number  of 
unknowns  is  generally  very  large.  Taking  uncertainties  into  account,  several  probabilistic 
approaches,  e.g.,  perturbation  method,  Kriging  predictor  and  random  matrix  theory,  have  been 
suggested  to  characterize  the  underlying  property  of  engineering  structures.  The  Bayesian 
inference  approach  has  shown  certain  advantages.  For  instance,  it  can  specify  the  model 
parameters  with  prior  information  in  the  form  of  probability  density  function  (PDF),  which  may 
be  viewed  as  imposing  soft  physical  constraints  to  enable  a  unique  and  stable  solution. 
Moreover,  this  approach  allows  the  computation  of  any  type  of  statistics  of  the  model  parameters 
to  be  identified. 


In  this  part  of  research,  we  develop  a  new  robust  algorithm  to  conduct  damage  localization 
and  identification  based  on  piezoelectric  impedance/admittance  information.  In  particular,  this 
algorithm  is  built  upon  the  Bayesian  inference  framework,  which  employs  forward  analysis- 
based  approach  instead  of  inversion-based  identification  procedures.  To  facilitate  such  analysis, 
the  tunable  inductive  circuitry  is  considered,  in  which  the  identification  performance  under 
different  inductance  tunings  is  evaluated  and  improved. 

When  applying  Bayesian  theorem  for  structural  model  updating,  the  hypothesis  0  is 
interpreted  as  the  vector  of  parameters  that  need  to  be  identified.  D  denotes  the  measured 
signature,  which  in  this  study  is  the  electrical  admittance  of  piezoelectric  transducer.  M  denotes 
modeling  assumptions,  reflecting  the  existing  experience  and  knowledge.  We  then  have 


M0|d,m)  = 


p{D\Q,M)p{^\M) 

|/7(D|0,M)/7(0|M)t/0 


(20) 
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The  prior  distribution  p{id  \  M)  expresses  the  initial  knowledge  of  eoncerned  parameters,  e.g., 
stiffness,  mass,  and  possible  damage  loeation.  The  choice  of  this  distribution  depends  on  how 
much  information  of  the  system  is  known.  In  this  research,  the  parameters  are  only  specified 
within  certain  range  based  on  the  prior  information.  Therefore,  for  the  sake  of  simplicity,  this 
term  can  be  defined  as  a  multivariate  uniform  distribution.  The  posterior  distribution  p{Q  \  D,M) 
indicates  the  updated  knowledge  of  the  parameters  0  conditional  on  the  prior  knowledge  and 
measured  admittance  information.  The  likelihood  function  /?(D|0,M)is  used  to  evaluate  the 
agreement  between  the  measurements  and  associated  model  output.  Specifically,  considering 
that  uncertainties  exist  in  real  measurement,  here  we  define  the  likelihood  function  as  a 
multivariate  normal  distribution  to  conduct  the  screen  of  model  output  over  0  space. 

.  (D-D(e)f  E-'(D-D(e)) 

=  ^  (21) 

V(2^)  |2| 


where  D  is  a  measured  admittance  vector  having  length  k  ,  and  D(0)  is  the  model  output 
parameterized  by  0  .  Z  is  the  covariance  matrix  of  D.  Under  this  framework,  the  damage  status 
of  the  structure  monitored  can  be  eventually  identified. 

Here  we  seek  to 
localize  and  identify 
small-size  damage  in  a 
plate  structure.  The  host 
structure  is  an  aluminum 
square  plate 

(  0.5m  X  0.5m  X  0.005m  ) 

attached  with 

piezoelectric  circuitry. 

The  plate  is  clamped  at 
its  two  ends,  and 
modeled  using  plate 
finite  elements.  The 
piezoelectric  transducer 

is  a  square  one  with  length  0.0025m,  which  covers  25  (5x5)  plate  elements.  An  inductor  is 
connected  to  the  transducer,  and  its  value  is  adjusted  to  3.6662  H,  in  order  to  generate  desired 
resonant  effects.  Assume  damage  occurs  in  a  single  element  with  40%  Young’s  modulus 
reduction,  and  is  located  near  the  center  of  the  plate.  As  no  prior  knowledge  of  possible  damage 
severity  (reduction  of  element  Young’s  modulus)  and  location  is  known,  we  choose  bivariate 
uniform  distribution  as  the  prior  distribution.  The  number  of  samples  parameterized  from  such  a 
distribution  is  equal  to  the  product  of  the  possible  damage  severity  candidates  and  damage 
location  candidates.  In  real  situation,  the  damage  severity  candidates  are  continuously  sampled 
within  certain  range,  and  damage  can  probably  occur  in  each  element  of  plate,  which  yields  a 
large  number  of  samples  based  on  prior  distribution.  Based  on  the  essence  of  Bayesian  inference, 
repeated  evaluation  of  the  likelihood  function  corresponding  to  such  samples  using  Monte  Carlo 
simulation  can  be  conducted,  where  evaluation  of  likelihood  function  depends  on  the  harmonic 
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response  analysis  with  multiple  exeitation  frequency  components.  For  the  purpose  of  illustration, 
here  we  consider  a  region  of  the  plate  where  damage  is  assumed  more  likely  to  occur  (including 
100  elements  marked  by  red  circle  in  Figure  10).  Moreover,  the  possible  damage  severity 
candidates  are  reduced  to  a  vector  [40%  60%  80%].  The  interested  frequencies  range  from 
1060Hz  to  1140Hz,  generating  uniformly  distributed  50  frequency  components.  Although  the 
damage  severity  does  not  continuously  vary,  we  can  use  meta-modeling  technique  such  as 
Gaussian  process  to  interpolate  and  predict  the  posterior  probabilities  under  unobserved  damage 
severity  candidates  in  the  future  study.  Taking  uncertainties  into  account,  the  0.5%  admittance 
measurement  standard  deviation  is  introduced  to  establish  the  likelihood  function.  Furthermore, 
in  this  case  study,  we  demonstrate  the  results  with  respect  to  the  different  inductance  values 
(z,  =  3.6662  and  =  3.5662  ),  aiming  at  facilitating  the  identification  process  and  further  providing 
the  adaptive  inductance  tuning  guideline  for  practical  experimental  manipulation. 

The  posterior  distributions  obtained  by  using  different  inductance  configurations  can  be  seen 
in  Figure  11.  As  such  distribution  inherently  is  high  dimensional,  we  convert  it  into  several  low 
dimensional  descriptions  in  order  to  highlight  the  comparison.  It  is  worth  mentioning  that  all 
the  identified  probability  distributions  in  demonstrated  figures  are  normalized.  Based  on  this 
series  of  results,  some  observations  can  be  summarized  as  follows. 
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(c)  Damage  index 


Figure  11  Posterior  distribution  versus  damage 
severity  and  damaged  eiement  index:  (a)  with 
I, ,  (b)  with ’  (^)  combination  of  Z, and  L^. 
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Damage  severity 

Figure  12  Posterior  distribution  of  damage 
iocation  under  40%  severity:  (a)  with  ,  (b) 

with £2,  (c)  combination  of  £jand  L^. 


Under  damage  severity  40%,  the  identified  distribution  by  using  the  information  of  both  Zj 
and  Z2  becomes  narrower,  which  indicates  better  identification  performance  as  compared  to  that 
obtained  by  using  admittance  of  individual  Z,  or  Z^  alone.  However,  under  damage  severity  60% 
or  80%,  this  identification  enhancement  is  not  that  obvious.  The  reason  is  that  the  nominal 
damage  severity  for  this  case  analysis  is  indeed  pre-specified  as  40%.  In  other  words,  other 
unreliable  damage  scenarios,  e.g.,  damage  severities  60%  or  80%,  cannot  yield  the  increasing 
probability  values  even  more  admittance  information  is  utilized.  This  is  certainly  in  agreement 
with  the  underlying  physical  characteristics.  The  advantages  of  using  the  concept  of  tunable 
inductance  circuitry  are  further  analyzed  by  comparing  the  probabilities  numerically.  When 
employing  the  tunable  inductance  circuitry,  several  high  fidelity  scenarios  have  the  increased 
likelihood  (larger  probability  values),  which  clearly  validates  the  feasibility  of  this  underlying 
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Probability 


idea.  The  identification  performance  under  Zj  is  better  than  that  of  ,  which  indicates  that  the 
inductance  tuning  technique  plays  an  important  role  in  conducting  damage  identification  process. 


Damage  severity  (c)  Damage  index 

Figure  13  Posterior  distribution  of  damage 
iocation  under  60%  severity:  (a)  with  1^ ,  (b) 
withZj’  (c)  combination  of  Xjand  Xj. 
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Damage  severity  (c)  Damage  index 

Figure  14  Posterior  distribution  of  damage 
iocation  under  80%  severity:  (a)  with  ,  (b) 
withXj,  (c)  combination  of  I^and  l^. 
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Interactions/T  ransitions 

This  research  is  directly  relevant  to  AEOSR’s  mission,  since  the  results  can  be  applied  to 
various  Air  Eorce  systems,  such  as  aircraft  fuselage,  space  structures,  and  gas  turbine  engines.  The 
Pis  have  had  various  interactions  and  technical  discussions  with  Dr.  Mark  Derriso  at  the  Wright- 
Patterson  Air  Eorce  Research  Eab  (AFRL)  on  issues  regarding  experimental  set  ups  and  DURIP 
ideas.  Dr.  Wang  also  has  had  communications  with  researchers  at  NASA  Glenn.  Dr.  Tang  has 
had  extensive  discussions  with  Pratt  &  Whitney  engineers  and  GE  Global  Research  Center 
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researchers  to  gather  engineering  insights  on  sensing  mechanisms  that  have  greatly  enhanced  the 
level  of  research. 
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Dr.  K.  W.  Wang  is  the  holder  of  the  Stephen  P.  Timoshenko  Collegiate  Chair  in  Mechanical 
Engineering  at  the  University  of  Michigan.  He  is  a  Fellow  of  the  ASME,  the  Institute  of  Physics, 
and  the  American  Association  for  the  Advancement  of  Science.  He  has  received  numerous 
recognitions  for  his  accomplishments;  such  as  the  SPIE  Smart  Structures  and  Materials  Eifetime 
Achievement  Award,  the  ASME  Adaptive  Structures  and  Materials  Systems  Prize,  the  ASME  N. 
O.  Myklestad  Award,  the  ASME  Adaptive  Structures  and  Material  Systems  Best  Paper  Award, 
the  ASME  Rudolf  Kalman  Best  Paper  Award,  the  NASA  Tech  Brief  Award,  and  the  SAE  Ralph 
Teetor  Award.  A  recent  paper  by  Jinki  Kim,  R.  E.  Harne,  and  K.  W.  Wang  (“Enhancing 
Damage  Identification  Robustness  to  Noise  and  Damping  using  Integrated  Bistable  and  Adaptive 
Piezoelectric  Circuitry,”  ASME  International  Design  Engineering  Technical  Conferences  & 
Computers  and  Information  in  Engineering  Conference,  DETC-20 14-34993,  Buffalo,  NY,  2014) 
has  won  the  Second  Place  Best  Student  Paper  Award  at  the  ASME  Conference  on  Mechanical 
Vibration  and  Noise. 
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